"0","richness <- rarefy(scpdSM, sample=minN, se=T) # because I specified *se=T*, this function outputs both S and se"
"0","S <- richness[1,]"
"0","S.se <- richness[2,]"
"0","# plot the result through time using base R plotting function"
"0","plot(S ~ ages, type=""l"", xlim=c(21000, 0), ylim=c(0, floor(max(S+S.se))+1), xlab=""Age (cal BP)"", ylab=expression('Richness' %+-% 'se'))"
